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are probably too faint to be a promising target for IACT observations, even with the prospective 
Cherenkov Telescope Array (CTA). 
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1. Introduction 

Revealing the nature of dark matter (DM) constitutes one of the most challenging tasks for 
astrophysics and cosmology. The most promising candidate to explain the numerous observational 
evidence for DM is a non standard model weakly interacting massive particle (WIMP), which self- 
annihilates into standard model particles [1]. As the WIMP mass is expected to be in the range of 
GeV to TeV, it is supposed to produce secondary y-rays in the same energy regime. This opens 
the window for an indirect detection of DM by y-ray telescopes. Various techniques exist for the 
detection of astrophysical y-rays. This contribution is focused on imaging atmospheric Cherenkov 
telescopes (IACT), in particular the prospective Cherenkov Telescope Array (CTA) [2], 

So far, no clear evidence for particle DM, neither direct nor indirect, has been found. Several 
targets have been considered for the indirect detection of DM, among them the Galactic center, 
neighboring dwarf spheroidal galaxies (dSphG), and galaxy clusters [3]. Extensive numerical sim¬ 
ulations in the framework of ACDM cosmology have shown that DM is also expected to form 
density structures on sub-galactic scales [4, 5]. Such subhalos in our Galaxy may constitute an 
interesting additional target for indirect DM searches. First, they are expected to lack of any con¬ 
taminating non-DM y-ray emission, resulting into a signal-to-noise ratio comparable or better than 
for the dSphG observations. Second, as the flux of secondary y-rays is expected to be proportional 
to the DM density squared, DM subhalos may cause a comparably high small-scale fluctuation of 
the isotropic y-ray background (IGRB), detectable and distinguishable from other sources [6]. 

Many studies have been performed to assess the y-ray fluxes from Galactic DM subhalos, 
by using the results of numerical simulations, e.g., [7, 8]. However, uncertainties about the spa¬ 
tial component of the signal, the astrophysical J-factor, still remain. The resolution of numeri¬ 
cal N-body simulations is limited to larger subhalo masses than relevant for a detectable signal. 
Therefore, extensive extrapolations and model assumptions have to be applied over many orders 
of magnitude. In this contribution, these model systematic uncertainties are studied with the new 
version of the semi-analytic code CLUMPY [9, 10]. Four different N-body simulation fit results and 
extrapolations are varied, and an uncertainty band is derived for the 7-factors and the fluctuations 
in the IGRB from DM subhalos in the Milky Way. 

2. y—rays from Galactic dark matter subhalos 

The y-ray flux S from annihilating dark matter at redshift z = 0 can be expressed as 



(2.1) 


where m x is the dark matter particle’s mass, (<rv) the velocity-averaged annihilation cross section, 
bj is the branching ratio for annihilation channel i, and (dNy/dEy)i the final-state y-ray spectrum 
in this channel. The second term, the 7-factor, 7 = // p 2 dl dD. is the integral of the squared DM 
density, p, over the line of sight, l, and AD = 2 jt( 1 — cos 6), the solid angle over the instrument’s 
aperture radius, 6. The flux S is proportional to 7, which is independent of any particle physics 
model. Therefore, in the following all quantities are expressed in terms of 7 instead of S. 

The mean density p of the total Galactic DM host halo is modeled by an Einasto profile, 
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on the basis of a best fit to the Aquarius simulations presented in [8], with the scale radius r s = 
15.14kpc and cxr = 0.17. The parameter I /(Xf is a measure of the steepness of the density profile 
close to r = 0, with r the distance from the halo center. Motivated by recent astronomical measure¬ 
ments [11], the halo is normalized by a local dark matter density of 0.4GeVcm~ 3 at r 0 = 8.0kpc. 
For a virial radius of 260kpc of the host halo, this results in a total halo mass of 1.1 • 1O 12 M 0 . This 
mass is in agreement with observations of stellar tracers [11] and is used in the following. 

The simulations indicate that the distribution of subhalos in the host halo follows an anti¬ 
biased distribution with respect to the unclustered, smooth DM distribution [4, 5]. This implies that 
most of the clustered mass is located in the outer parts of the host halo. For the model variation, 
two parametrizations of the substructure distribution ~p sub within the host halo are compared: The 
Aquarius-simulation [4] suggests an Einasto profile for the subclustered mass with (Xf = 0.678 
and r s = 199kpc (model EINASTO). The Via-Lactea II subhalos [5] are described by a slightly 
less anti-biased distribution, modeled by the number distribution N su b{< r) given in [12] (model 
MADAU). All subhalos are described by an Einasto profile with the same (Xu as for the total halo. 

The presence of subhalos with masses in the range [1O~ 6 M 0 , 1O 1O M 0 ] is assumed [13] 1 . The 
simulations [4, 5] indicate a power-law function dN su b /d M su b ~ M~ b "' of the subhalo distribution, 
with an index 1.9 < a m < 2.0 for the resolved subhalos with masses M sub > 10 5 M & . The extrap¬ 
olation of dN su b/dM su b down to the lowest subhalo masses is very sensitive to (X m . Therefore, the 
results for the two indices a m = 1.9 and a m = 2.0 are compared. A fixed number of 150 massive 
dSphG-like subhalos in the mass range [1O 8 M 0 ,1O 1() M 0 ] is supposed, in agreement with [4]. 

An important property of a subhalo is the concentration c(M su b ) = R\/(s,sub, which connects 
its scale radius r S)SU b with its virial size R\ as a function of its mass M su b • As the annihilation 
y-ray emission is proportional to the density squared, a higher concentration significantly boosts 
the brightness of a subhalo. The rather optimistic (highly concentrated) model from [14] (model 
BULLOCK) is compared to the conservative one (flatter subhalo concentrations) from [13] (model 
SANCHEZ). Finally, the authors of [14] find a significant log-norm scattering a c = A(logc) around 
the mean concentrations c{M su b), which also results into a net boost of the signal. Their result for 
hosted subhalos, o c = 0.24, is compared to a more conservative assumption of a, = 0.14. 

From the two values for each of the four varied quantities, 16 benchmark models are created, 
listed in Tab. 1 . Model 9 is chosen to be the reference model with an average flux yield. 
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Ctc 

0.14 

0.14 

0.14 

0.14 

0.14 

0.14 

0.14 

0.14 

0.24 

0.24 

0.24 

0.24 

0.24 

0.24 

0.24 

0.24 

aj 

2.05 

2.06 

2.15 

2.15 

2.03 

2.02 

2.12 

2.13 

2.04 

2.05 

2.14 

2.15 

2.03 

2.04 

2.09 

2.12 


Table 1: Parameters for the 16 models investigated in this study. The reference model is marked in 
gray. S/B stands for the SANCHEZ/BULLOCK c(M su b ) parametrization, E/M: p mb parametrization by 
an EINASTO/MADAU profile. Bottom row: Fit values for the exponent ay of a power law fit for 
dN su b/dJ ~ J aj in the regime J < 10 17 GeV 2 cm . The fitting uncertainty is o aj < 5 ■ 10 2 for all models. 


1 Subhalos are drawn by CLUMPY down to a user-defined accuracy in terms of the 7-factor contribution from subha¬ 
los compared to the background emission from the unclustered “smooth” DM in the halo. For all models, all subhalos 
yielding J > 10 16 GeV 2 cm~ 5 are correctly simulated in this study. 
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3. Setup of the simulations 

The Galactic DM halo is simulated for a circular field of view (FOV) centered at the Galactic 
pole b = 90° with diameter d = 120°. This FOV covers a fraction /^ y = 25% of the sky. The choice 
of the FOV is motivated by the projected extragalactic sky survey achievable with CTA [15], an 
area where it can be searched also for Galactic DM subhalos. To take into account the statistical 
variation of the observed sky in terms of cosmic variance, each simulation is rerun 10.000 times. 
In the same way, the search for DM subhalos within a single FOV of d = 10° is investigated, using 
a separate simulation. 

4. Results 

In Fig. 2, the average integrated source count distribution is shown, defined as 

co 

N sub {>J) = (4.1) 

j 

Here, N su b(> J ) represents the number of subhalos on the chosen patch of the sky brighter than a 
given 7-factor. The DM subhalo source count distribution dN sl ,b /d7 can be described by a power- 
law according to dN su b/dJ ~ 7 aj in the regime 7 < 10 17 GeV 2 cm 5 . In Tab. 1, the fit results 
of the exponent otj are listed for each of the 16 models. These results indicate that the dN su b/dJ 
distribution is dominantly determined by the mass distribution d N su b /d M su b M- U b m - 

The integration angle, within which the 7-factor of each subhalo is calculated, is set to 9\ = 
0.1° to illustrate the power-law behavior of the source count distribution. In Fig. 1, it is shown 
for model 9 that the largest part of most of the subhalos is encompassed within 0\. Similar 
results are obtained for all models. As the annihilation y-ray emission is strongly peaked at 
the subhalo center, only the signal of the most extended subhalos gets significantly truncated 
by ©i. Fig. 1 suggests a correlation between the apparent size and the brightness of the sub¬ 
halos. Therefore, a cut-off of the power-law distribution at high 7 is visible in Fig. 2. Contin¬ 
uous emission from the unclustered DM (“smooth halo”) and other background decreases pro¬ 
portional to 6 2 , faster than the subhalos’ flux [16]. Consequently, a smaller signal region might 
improve the signal-to-noise ratio in a real experiment. 

Only a slight correlation between the subhalo 
masses and their 7-factors has been found for all 
models (not shown). From this follows that distant 
massive dSphG-like subhalos and nearby low-mass 
subhalos similarly account for the highest 7-factors. 

For a given DM annihilation channel, cross sec¬ 
tion and particle mass, the fluxes corresponding to the 
7-factors are plotted in Fig. 2. For the annihilation 
channel XX bb, CTA will be most sensitive to a DM 
particle mass of 500 GeV < m x < 1 TeV [17]. For this 
mass range, the annihilation cross section is already 
constrained to (ctv) < 3 • 10~ 25 cm 3 s _1 by observa¬ 
tions of the dSphGs by the Fermi-LAY instrument [18]. 



1(T 2 KT 1 1 

a s =atan(r s /D) [deg] 


Figure 1: Relation between the angular 
size a s and the brightness J(9\ = 0.1°) of 
the subhalos for the survey region & model 
9. D is the subhalo distance from the ob¬ 
server. The dashed line indicates a s = 6\. 
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flux > 100 GeV in Crab units for x —> bb , ( av ) = 3-10 25 cm 3 s m x = 500 GeV 
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Astrophysical J-factor in /ogio(J/GeV 2 cm 5 ) 


Figure 2: Integrated source count distribution, N su b, for all 16 benchmark models on f s b y = 0.25 and for 
01 = 0.1°. The dependence of /V v „/, on the /-factor is shown on the lower abscissa. On the upper abscissa, 
the dependence of N su b on the flux is shown, for the assumed particle physics model. The continuous/dashed 
lines distinguish the SANCHEZ/BULLOCK concentration model. The corresponding /-factor from Sculp¬ 
tor dSphG for 0i = 0.1° [19] is plotted for comparison. 



Zo3i 0 (//GeV 2 cm- 5 ) 

Figure 3: N st ,b of the four models that differ only 
by the depicted parameter from model 9 (cr t = 0.24, 
p sub MADAU, a m = 1.9, c(M suh ) SANCHEZ), nor¬ 
malized to N su b °f model 9. The dashed line 
marks the / of the brightest subhalo expected 
(N su b(> /) = 1) for model 9. 



Zopio(//GeV 2 cm s ) 


Figure 4: N su b f° r model 9 and a FOV diameter of 
d = 10° at different distances A b from the Galactic 
center. The maximum difference between all solid 
curves is < 10%. The dashed lines denote the mean 
emission from the smooth galactic DM component 
for each FOV, and 0i = 0.1°. 


Thus, in the range of the considered models and a coverage of f s ty = 25%, a survey sensitivity to 
some 10 4 of the Crab Nebula flux is necessary to discover at least one subhalo above 100 GeV 
with IACTs, provided that (av) ~ 3 • 10~ 25 cm 3 s _1 . 

Fig. 3 shows the impact of the varied parameters on N su b compared to the reference model. For 
the brightest subhalos, the model uncertainty is generally < 50%. An exception from this occurs 
when switching from the SANCHEZ to the BULLOCK concentration model for c(M su b). In this 
case, the number of brightest subhalos is increased by more than a factor of three. 
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Fig. 4 shows the simulation results for N su b within a single CTA-like FOV of d = 10°. Here, 
N su b is given for different distances A b from the Galactic center, perpendicular to the Galactic plane. 
Two conclusions can be drawn from this. First, Fig. 4 indicates that at r Q , the subhalos appear 
almost isotropically in the sky. Consequently, the relation N su b(AQi) /N su b(AQ. 2 ) ~ AG|/AFL 
holds for arbitrary positions in the sky (and is confirmed by comparing Fig. 2 and 4). However, 
pointing close to the Galactic center heavily decreases the contrast between the subhalos and the 
smooth DM halo emission, the latter indicated by the dashed lines. Second, a sensitivity to 7-factors 
of 7 < 10 18 GcV 2 cm 5 is needed to resolve at least one subhalo within a given single FOV. Such 
7-factors are about a factor 10 smaller than those of the most promising dSphGs, where observation 
times of A> 500hrs are expected to reach competitive limits ( ov ) ~ 10 _25 cm 3 s _1 with CTA [17]. 

In Fig. 5, the fluctuation angular power spectra (fluctuation APS) of the subhalos are shown, 
comparing two of the four varied parameters against the reference model. The fluctuation APS is 
a measure of the relative fluctuation, or small scale anisotropy, of a signal, which can be possibly 
detected above a much larger background holding less fluctuation. The fluctuation APS C F of a 
map /(#, (p), defined on a patch f^y of the sphere (7 = 0 elsewhere) is given by 


r F 


1 

fsky 7 2 (27 + 1) 


+£ 

L i 12 

| dim | 5 

m——£ 


(4.2) 


with the spherical coordinates $ and <p, and I the mean value of / on the defined domain. The a\ m 
are the coefficients of the map decomposed into spherical harmonics Y\ m , 

t-max 

I($,(p) = Y Yj a hn Y lm ($,(p). (4.3) 

£—0 m——£ 


For Fig. 5, an integration angle 6b = 0.015° is chosen, which corresponds to half the value of 
the m a ximum resolution projected for CTA at energies above 1 TeV [15]. 

For a FOV with diameter d < 10° and an angular resolution of CipsF > 0.03° [15], CTA could 
be sensitive to multipoles 100 < 7 < 1000 [6]. As shown in Fig. 5, in this regime, the fluctuation 
APS remains largely constant for all models, with only a mild lowering of the amplitudes Cf 
towards higher multipoles. Only in the case of the variation of a m from 1.9 to 2.0, a drastic change 
for the fluctuation APS results. This is shown in Fig. 5 (left). For a m = 2.0, a substantially higher 
amount of total mass is bound into subhalos. Although this results into a boost of the overall 
signal, the effective fluctuation heavily decreases because of the large background of the low-mass 
subhalos. For all other cases, the fluctuations of the subhalos remain largely the same. This applies 
even for the significant signal boost for the BULLOCK case, which is shown in Fig. 5 (right). 


5. Conclusions 

Different models of the abundance and properties of Galactic DM subhalos have been com¬ 
pared. Within the varied parameters, the source count distribution dN su b/dJ of these subhalos has 
been constrained to an uncertainty band of one order of magnitude. The modeling of the inner 
concentration c(M su b) causes the biggest systematic uncertainty in the number of detectable subha¬ 
los. The subhalos appear nearly isotropically distributed in the sky, allowing a telescope to point as 
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Figure 5: Fluctuation APS Cf of the subhalos on f\ky = 0.25, comparing the variation of the index a m (left) 
and the c(M su b) model (right) against model 9. The fluctuation APS of the smooth halo on the simulated 
patch of the sky is plotted for comparison. The bands indicate the la containment ranges around the mean 
spectrum of the simulated galactic halo, assuming a log-normal variation of the Cf. Note that for the left 
panel, also the smooth contribution varies between the compared models. However, this difference in C[ 
between the smooth contributions is smaller than the resolution of the plots and not distinguished here. 

far as possible away from the Galactic plane to reduce the background from the smooth DM halo 
emission and diffuse astrophysical foregrounds. 

However, for all benchmark models, Galactic subhalos are probably too faint to be a promising 
target for current and next-generation IACTs, provided the projected available observation times. 
Still for a sensitivity to a fraction of 10 4 of the Crab Nebula flux above lOOGeV, a detection 
of individual subhalos might be possible only for cross sections (<rv) > 10 -25 cm 3 s -1 , which are 
already in tension with recent measurements. Pointed observations of the Galactic center and the 
dSphGs with CTA promise competitive limits on the annihilation cross section, probably requiring 
much less observation time. 

At the same time, it has been shown that a source count distribution with aj < 2.1 results 
in a fluctuation APS up to the order of 10 3 < Cf < 10 2 in the multipole range of interest of 
Cherenkov telescopes. In principle, such a comparably high fluctuation might distinguish dark 
matter from other, more dominant unresolved source classes. Possibly, even more competitive lim¬ 
its on the annihilation cross section can be obtained than by the search for individual subhalos [6]. 
However, given the low intensity of the overall signal from Galactic DM subhalos, these limits are 
most likely still not competitive with those obtainable from (combined) pointed observations. 

We did not assess the extragalactic contribution to DM induced fluctuations of the IGRB. 
However, the results of [8] indicate that the extragalactic emission is expected to be similar to 
the Galactic signal in magnitude and in angular scales, suggesting a similar sensitivity to IGRB 
fluctuations as for the Galactic subhalo case. Still, it should be emphasized that all indirect DM 
targets suffer from major systematic uncertainties. With this in mind, the search for fluctuations in 
the IGRB should not be neglected as a complementary strategy to unveil the nature of DM. 
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